Abstract

Violence against women, especially in the form of sexual abuse and assault, is a serious issue for public health and an infringement on women’s human rights. In India, Uttar Pradesh has the highest incidence of all crimes against women; rape and dowry deaths are two of the most common crimes against women. The study aimed to gain insights into the spatial and temporal dynamics of these crimes to pinpoint effective countermeasures using Bayesian models. Data from the NCRB was used to analyze the prevalence of rape and dowry deaths across the 70 districts of Uttar Pradesh from 2001 to 2014. The hierarchical structure of the data and its spatial and temporal characteristics were taken into account using spatial and temporal modeling techniques. The results showed the highest incidence rate of the two crimes in the central and western parts of the state, with an overall decrease from 2001 to 2003, followed by an increase in the following years. The study emphasized the need for effective policy and intervention measures to reduce and prevent violence against women in Uttar Pradesh, focusing on the areas where the incidence of these crimes is highest.

Introduction

Sexual assault and romantic partner violence are serious public health issues and human rights violations. According to the World Health Organization (WHO), around one in three (30%) women worldwide have experienced physical, sexual, or both forms of intimate relationship abuse or non-partner sexual violence at some point (WHO 2021). The National Crime Records Bureau (NCRB) stated that Uttar Pradesh, the most populous state in India, has the highest rate of overall crimes against women in 2017. In 2002, the Protection from Domestic Violence Bill was established to “protect the rights of women who are victims of any kind of violence in the family and to provide for things that are related to or related to that violence (Sabha 2002).” Rape incidences and dowry deaths are the two most common crimes against women in Uttar Pradesh. These crimes are distressing and harmful to women’s physical and emotional health.

To learn more about the spatial and temporal dynamics of rape incidences and dowry deaths, we conducted an analysis of the prevalence of crimes against women in Uttar Pradesh. The analysis emphasizes the importance of considering spatial and temporal variables when investigating crimes against women, laying the groundwork for future studies on the underlying social, economic, and cultural factors that influence the prevalence of these crimes in Uttar Pradesh and other regions of India. The purpose of this analysis is to learn more about how these incidents have changed over time and across different districts in the state. The main goal of the analysis is to help guide policy and intervention efforts that aim to reduce and stop violence against women.

Methods

Data

The data was obtained from the NCRB, and spatial and temporal modeling techniques were used to analyze patterns and trends in the incidence of rape and dowry deaths in Uttar Pradesh. According to the NCRB, the data was collected by the State Crime Records Bureaux (SCRBx) from the District Crime Records Bureaux (DCRBx) and sent to the NCRB at the end of each calendar year (Bureau 2022). The data contained information on rape incidences, dowry deaths, the female population between 15 and 49 years old during the period 2001–2014, the standardized incidence/mortality ratio for rapes and dowry deaths, and the IDs for the area, year, and area-time identifiers. The total female population between 15 and 49 years old in 2001 was 36,343,890, and in 2014 it was 50,563,541. Fig. 1 shows a map of Uttar Pradesh with the outline of the districts, and Table 1 in the supplementary material shows the ID for the districts to help identify areas on the map. Figure 2 displays the trend for the crude rates over the years for Uttar Pradesh, which dropped sharply from 2001 to 2004 but rose in the end for both rape incidences and dowry deaths.

Table 2 shows the descriptive statistics for the rape incidence and dowry death cases in Uttar Pradesh, India, over the years. From 2012 to 2014, the mean for rape incidence increased dramatically. The dowry death incidences rose from 2004 to 2009, then stabilized until 2014. Both tables showed decreasing trends from 2001 to 2003 and 2005.

Table 2: Descriptive statistics
year Minimum_R Q1_R Mean_R sd_R Q3_R Maximum_R Minimum_D Q1_D Mean_D sd_D Q3_D Maximum_D
2001 1 (Sant Ravidas Nagar Bhadohi) 13.00 27.94 21.52 41.00 93 (Bareilly) 4 (Shrawasti) 18.00 31.57 19.00 43.75 88 (Kheri)
2002 0 (Sant Ravidas Nagar Bhadohi) 9.00 20.19 14.73 30.75 73 (Sitapur) 3 (Shrawasti) 14.25 27.04 18.08 34.75 83 (Sitapur)
2003 0 (Sant Kabir Nagar) 5.00 12.97 11.06 19.50 47 (Sitapur) 3 (Balrampur) 10.25 18.89 11.49 24.00 55 (Agra)
2004 3 (Mau) 9.25 19.93 14.98 25.75 72 (Sitapur) 3 (Balrampur) 14.25 24.40 15.34 29.00 71 (Sitapur)
2005 1 (Lalitpur) 7.00 17.34 14.17 24.00 61 (Aligarh) 1 (Lalitpur) 12.25 22.33 13.93 26.75 70 (Sitapur)
2006 2 (Sant Ravidas Nagar Bhadohi) 9.00 18.76 12.14 26.00 51 (Lucknow) 7 (Sant Kabir Nagar) 14.25 25.69 14.37 34.75 67 (Kanpur Nagar)
2007 1 (Sant Ravidas Nagar Bhadohi) 10.00 23.51 16.57 32.50 82 (Lucknow) 4 (Shrawasti) 16.00 29.63 17.28 36.75 78 (Agra)
2008 2 (Sant Ravidas Nagar Bhadohi) 12.00 26.70 19.03 35.75 82 (Lucknow) 5 (Balrampur) 17.25 31.96 18.73 38.75 88 (Aligarh)
2009 3 (Chandauli) 13.00 25.11 17.52 35.25 77 (Bareilly) 8 (Lalitpur) 19.25 31.87 17.98 40.75 83 (Agra)
2010 1 (Sant Ravidas Nagar Bhadohi) 10.25 21.89 17.37 26.00 75 (Lucknow) 5 (Balrampur) 18.25 31.44 19.67 40.00 95 (Aligarh)
2011 2 (Sant Ravidas Nagar Bhadohi) 14.25 29.14 20.63 39.00 89 (Moradabad) 6 (Sant Kabir Nagar) 17.00 33.16 18.69 41.75 95 (Aligarh)
2012 4 (Mirzapur) 15.00 28.00 17.39 35.75 86 (Bareilly) 5 (Balrampur) 19.00 32.03 17.88 40.75 97 (Aligarh)
2013 5 (Mirzapur) 23.25 43.53 28.52 53.75 119 (Ghaziabad) 5 (Sant Kabir Nagar) 19.00 33.31 19.52 41.25 98 (Agra)
2014 5 (Shrawasti) 23.00 49.46 31.74 69.00 164 (Ghaziabad) 6 (Sonbhadra) 23.25 35.26 18.39 46.75 98 (Aligarh)

Spatial model

The hierarchical Poisson log-linear model was used to estimate the incidence rates of rape and dowry deaths separately while accounting for the hierarchical structure of the data and controlling for spatial and temporal dependencies in the data. Specifically, the model was used to estimate the relative risks of rape incidence and dowry deaths across the districts of Uttar Pradesh, taking into account the spatial correlation between neighboring districts. The BYM2 prior was employed to model the spatial random effects, which combines a conditional autoregressive (CAR) prior and an independent and identically distributed (iid) prior. This approach allowed for modeling both the spatially structured and unstructured random effects.The model is specified as: \[ \begin{eqnarray} O_i &\sim & \text{Poisson}(\rho_i E_i)\\ \eta_i &= & \log \rho_i = b_0 + b_i\\ \boldsymbol{b} &= & \frac{1}{\sqrt{\tau_b}}(\sqrt{1-\phi}\boldsymbol{v}_{*} + \sqrt{\phi}\boldsymbol{u}_{*})\\ \end{eqnarray} \] where \(\boldsymbol{v}_{*}\) and \(\boldsymbol{u}_{*}\) are standardised versions of \(\boldsymbol{u}\) and \(\boldsymbol{v}\). \(u_i\) is the spatially structured component defined by an intrinsic CAR prior: \(\boldsymbol{u}\sim ICAR(\boldsymbol{W}, \sigma^2_u)\), \(v_i\) the unstructured component defined with prior: \(v_s \overset{iid}{\sim} \text{Normal}(0,\sigma^2_v)\).

Spatio-temporal analysis

A spatio-temporal analysis using a separable space-time model with type 1 interaction was conducted to examine both the spatial and temporal patterns of the incidence of these types of violence. The incidence of rape and dowry may vary both spatially and temporally, understanding these patterns will allow us to identify areas and periods of high risk and target interventions effectively. A spatial-temporal analysis with a separable space-time model with type 1 interaction was used to identify areas and periods of high risk, as well as to test hypotheses about the factors that are driving the incidence. The model can be utilized to test whether there are spatial or temporal clusters of high incidence rates or whether certain socio-economic or demographic factors are associated with higher incidence rates.

\[ \begin{eqnarray} O_{it} &\sim & \text{Poisson}(\rho_{it} E_{it}) \\ \log \rho_{it} &= & b_0 + b_i + \gamma_t + \psi_t + \delta_{it} \\ \boldsymbol{b} &= & \frac{1}{\sqrt{\tau_b}}(\sqrt{1-\phi}\boldsymbol{v}_{*} + \sqrt{\phi}\boldsymbol{u}_{*})\\ \gamma_t & \sim & \hbox{RW(1)}\\ \psi_t & \sim & N(0,\sigma^2_{\psi})\\ \delta_{it} & \sim & \hbox{Normal}(0, \sigma^2_{\delta}) \end{eqnarray} \]

WAIC

WAIC was calculated for both the spatial and spatio-temporal models to compare their goodness of fit. It is used for model selection and comparison in Bayesian data analysis. The idea is to estimate the out-of-sample predictive accuracy of a given model using information from the posterior distribution. The spatio-temporal model yielded a better score, indicating that it provided a better explanation for the observed data. For both rape incidence and dowry deaths, the WAIC was about 10 times better for the spatio-temporal type 1 interaction model compared to the spatial model (WAIC for Rape: spatial-649, SpatTemp type1-6245; WAIC for dowry: spatial-659, SpatTemp type1-6335).

Results

The average SMRs for rape incidence and dowry deaths cases over the period 2001 - 2014 were mapped (fig 3 and fig 4). According to the figures, the grid (78,27) area has the highest SMRs compared to most of the other areas. The SMRs were relatively low for the west regions for both rape incidences and dowry deaths.

From the maps of the posterior mean of the residual relative risks and posterior probabilities, the posterior mean of the residual RRs map indicates that districts in the north and northeast have a relatively high risk of rape incidences, while the districts mainly in the east have a relatively low risk. For the dowry death cases, the posterior mean of the residual RRs map shows a relatively high risk for the north-central districts, while the east also has a relatively low risk in this case. The maps for the posterior probabilities for both rape incidence and dowry death cases suggest an association with some type of covariate that results a higher incidence of rape or dowry deaths. This indicates that the cases were not due to chance. However, we do not have enough information from the dataset on the covariates to specify which factor is driving the incidences.

Fig. 4a: RR and PP Spatial model for rape incidence cases

Fig. 4a: RR and PP Spatial model for rape incidence cases

Fig. 4b: RR and PP Spatial model for dowry death cases

Fig. 4b: RR and PP Spatial model for dowry death cases

The posterior summary of the hyperparameters for both rape incidence cases and dowry death cases both show that 38% of the spatial variability is explained by the spatially structured component, meaning there is some underlying spatial pattern in the data that can be accounted for by the model. This suggests that there are some spatial factors that are associated with the incidence of rape and dowry deaths in Uttar Pradesh.The spatially structured component of the model captures this underlying spatial pattern by modeling the spatial autocorrelation in the data. It assumes that observations that are closer together in space are more similar to each other than observations that are farther apart. However, the 38% of the spatial variability explained by the spatially structured component does not mean that spatial is the only important factor affecting the incidence of rape and dowry deaths. Other factors such demographic, socio-economic, or cultural factors, also influence the incidence of these types of violence.

Figure 5 presents the overall spatio-temporal trend, a huge decrease in risk is noticed for both rape incidences and dowry death cases from 2001 to 2003. However, the spatio-temporal risk increases, with the increase for dowry death cases stabilizing around 2010. This trend could be due to the establishment of a law against domestic violence in 2002.

Fig. 5: ST model Int I for Rape Incidences and Dowry Deaths

Fig. 5: ST model Int I for Rape Incidences and Dowry Deaths

Comparing the spatio-temporal model maps of the residual RRs and posterior probabilities for rape cases and dowry incidence with just the spatial models (fig. 6 and 7), the spatio-temporal model for rape incidence did not show difference in the spatial residuals, but the spatio-temperal model for dowry deaths showed a decrease in the spatial residuals.

Fig. 6: Spatio-temporal model: Map of the residual RRs and posterior probabilities for rape cases

Fig. 6: Spatio-temporal model: Map of the residual RRs and posterior probabilities for rape cases

Fig. 7: Spatio-temporal model: Map of the residual RRs and posterior probabilities for dowry death cases

Fig. 7: Spatio-temporal model: Map of the residual RRs and posterior probabilities for dowry death cases

Figure 8 illustrates that from 2001 to 2003, the maps for both rape incidence and dowry death become lighter suggesting a decrease in risk. However, from 2003 to 2008 the maps gradually darkened, and the geographical pattern appears to stabilize from 2010 onwards. Most of the low-risk districts are located in the eastern part of Uttar Pradesh, with a few additional low-risk districts in the northwestern corner. High-risk districts are primarily found in the western-central region of the state.

Discussion

The analysis of the prevalence of crimes against women in Uttar Pradesh, India, reveals the relevance of considering spatial and temporal variables. The results show that the highest rates are observed in the central and western part of the state. The temporal trend shows that the two crimes decreased from 2001 to 2003, then increased in the following years. The results suggests some spatial factors contribute to the state’s rape and dowry deaths, but it’s important to note that the spatially structured component of the model does not account for all of that are contributing to the incidence of these crimes, which need to be further investigated.

There were some limitations in the study. Fewer than 40% of the women who experience domestic violence speak up. Fewer than 10% of those who did ask for help also went to the police (WHO 2022). Estimates of unreported rape cases in India vary greatly. According to the National Crime Records Bureau’s 2006 report, approximately 71% of rape incidents go unreported (Indian Political Science Association 2009). The dataset did not include women’s age or age groups, making it difficult to identify age-specific patterns. Additionally, the dataset did not show any covariates, such as socioeconomic and demographic factors. Although the standardized mortality ratios might have been adjusted for covariates, we cannot be certain, which makes it difficult to draw definitive conclusions from the posterior probabilities. The analysis underscores the necessity for targeted policy and intervention efforts to mitigate violence against women in Uttar Pradesh. It emphasizes considering temporal trends when designing interventions. While the 2002 Protection from Domestic Violence Bill was a positive step, addressing the underlying social, economic, and cultural factors contributing to such violence remains crucial.

Violence against women is a pervasive issue with far-reaching consequences for individual women, their families, and society as a whole. The public health sector should take a multifaceted approach, which includes, self-defense classes to help women defend themselves physically but also fosters mental resilience and self-esteem. Community education to help change societal attitudes and promote a culture of respect and equality. Legal reforms to advocate for and supporting legal reforms can help to strengthen laws that protect women from violence.Providing therapy and counseling services to survivors of violence can help them heal from the emotional and psychological effects of their experiences.Offer comprehensive support services, such as shelters, hotlines, and financial assistance, can help women escape abusive situations and rebuild their lives. By implementing a combination of these strategies, the public health sector can contribute to reducing the prevalence of violence against women and improve the overall health and well-being of those affected by it.

Supplementary material

References

Bureau, National Crime Records. 2022. “Crime in India.” https://ncrb.gov.in/en/crime-india.
Indian Political Science Association. 2009. The Indian Journal of Political Science. Vol. 70.
Sabha, Rajya. 2002. “Hundred-Twenty-Fourth Report on the Protection from Domestic Violence Bill, 2002.” New Delhi, India: Parliament of India. http://164.100.47.5/rs/book2/reports/HRD/Report124th.htm#:~ :text=The%20Protection%20from%20Domestic%20Violence%20Bill%2C%20200 2%2C%20seeks%20to%20protect, connected%20therewith%20or%20incidental%20thereto.Rape.
WHO. 2021. “Devastatingly Pervasive: 1 in 3 Women Globally Experience Violence.” Geneva, New York: World Health Organization. https://www.who.int/news/item/09-03-2021-devastatingly-pervasive-1-in-3-women-globally-experience-violence.
WHO. 2022. “Facts and Figures: Ending Violence Against Women.” Geneva, New York: World Health Organization. https://www.unwomen.org/en/what-we-do/ending-violence-against-women/facts-and-figures.